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METHOD AND DEVICE AS WELL AS COMPUTER PROGRAM FOR CONTROLLING 

AN INTERNAL COMBUSTION ENGINE 

Background Information 

The present invention is directed to a method and a device as 
well as a computer program for controlling an internal 
5 combustion engine. 

To control an internal combustion engine, it is known from DE 
42 3 9 711 Al (U.S- Patent 5,558,178) to convert a setpoint 
value for a torque of the internal combustion engine into an 

10 actuating variable for influencing the air supply to the 
combustion engine, to set the ignition angle and/or to 
suppress or connect the fuel supply to individual cylinders of 
the internal combustion engine. In addition, the controlling 
of the fuel/air ratio to realize the predefined torque value 

15 is also known from WO-A 95/24550 (U.S. Patent 5,692,471). 

Furthermore, in the known approaches the actual torque of the 
internal combustion engine is calculated taking the actual 
engine settings into account (charge, fuel metering and 
ignition angle) . Among others, the engine speed, load (air 

20 mass, pressure etc.) and possibly the exhaust-gas composition 
are utilized. 

In the context of these calculations, a torque model for the 
internal combustion engine is used, which is utilized both for 
25 determining the actuating variables and also for ascertaining 
the actual values. The essence of this model is that values 
for an optimum torque of the internal combustion engine and 
for an optimum ignition angle are determined as a function of 
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operating points, which are then corrected via efficiency 
values in accordance with the actual setting of the internal 
combustion engine. 

5 To optimize this model, it is known from DE 195 45 221 Al 

(U.S. Patent 5,83 2,8 97) to correct the value for the optimum 
ignition angle as a function of variables influencing the 
efficiency of the internal combustion engine such as 
exhaust-gas recirculation rate, engine temperature, intake-air 
10 temperature, valve -overlap angle, etc. 

In practice, however, it has been discovered that this known 
approach may yet be optimized further, especially with respect 
to simplifying the application, optimizing the computing time 

15 and/or considering the working-point dependency of the 

correction of the optimum ignition angle, especially the 
dependency on the inert-gas rate. In some operating states, in 
particular, the known torque model does not deliver 
satisfactory results. Such operating states are especially 

20 states with high inert-gas rates, i.e., states having a large 
proportion of inert gas (by external or internal exhaust -gas 
recirculation) , which are caused by an overlap of intake and 
discharge-valve opening times and, above all, [occur] with 
small to medium fresh-air charges. Furthermore, these are 

25 operating states having high charge movements. Since these 

effects are not fully taken into account, the calculated basic 
variables make it impossible for the known procedure to obtain 
a precise torque calculation. 

3 0 Another problem is that the known model, while taking the 

essential influences on the ignition angle into account when 
determining the maximum torque, does not consider the 
different slopes of the correlation between ignition angle and 
the maximum engine torque which is able to be attained with 

35 different mixtures in the instantaneous operating point of the 
engine. In optimizing the known model, a compromise must thus 
be found, which usually consists in the optimum ignition angle 
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no longer coinciding with the ignition angle of the optimum 
torque. For example, in mixtures having good combustibility, 
which have a highly retarded ignition angle when generating 
the best torque, the optimum ignition angle is markedly more 
5 retarded than this ignition angle, whereas, in mixtures having 
poor ignitibility, it is considerably more advanced than the 
ignition angle at which the maximum torque is generated. This 
characteristic is sketched in Figure 3 . The dashed line there 
represents the curve provided by optimizing the known model. 
10 As can be seen, the actual and the model curve no longer 
correspond precisely. In one example of application, a 
resulting torque fault of up to 2 0% relative to the optimum 
torque MIopt has come about at the best possible optimum 
ignition angle. 

15 

Thus, it is obvious that the greater the variation in the 
steepness of the ignition hook (correlation between optimum 
ignition angle ZW^pt and torque MI with respect to top dead 
center of ignition phase ignition-TDC for an operating point 

20 of the engine) , the greater the deviation between the optimum 
ignition angle and the maximum ignition angle at which maximum 
torque is attained in mixtures that combust extremely well and 
those that combust poorly, and the greater the errors of the 
torque model. The known optimizer of the model parameters is 

25 unable to find a population of the model, in particular the 
optimum ignition angle, so that the torque model has low 
tolerance across the entire ignition- angle range. 

As a result, the torque model must be optimized further, in 
3 0 particular with respect to engine -control systems having high 
inert-gas rates, such as engine -control systems having 
variable valve timing and/or charge -movement flap. 

From the not pre-published German patent application 101 49 
35 477.7, a method and a device as well as a computer program for 
controlling an internal combustion engine are known, a torque 
model being used within the framework of calculating 
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instantaneous variables and/or actuating variables* In doing 
so, a basic value ascertained under standard conditions is 
corrected as a function of the inert-gas rate and/or the 
valve-overlap angle. Moreover, to further improve the 
5 precision of the model, the efficiency of the conversion of 
the chemical into mechanical energy by which the optimum 
torque value is corrected is determined as a function of the 
deviation between an optimum ignition angle and an 
instantaneous ignition angle as well as an additional variable 
10 that represents the combustion performance of the mixture, the 
latter being the optimum ignition angle in this case. 

Summary of the Invention 

15 The method according to the present invention, the device 
according to the present invention, the computer program 
according to the present invention and the computer-program 
product according to the present invention, having the 
features of the independent claims, have the advantage over 

2 0 the related art that the efficiency of the conversion of the 

chemical energy into mechanical energy is determined at least 
as a function of a variable characterizing the combustion 
center point and a variable characterizing the opening instant 
of a discharge -side gas-exchange valve. In this way, the 
25 reduction of the indicated engine torque associated with an 

advanced opening of the discharge- side gas-exchange valve, is 
taken into consideration in the torque model. As a result, it 
is possible to achieve high accuracy in the indicated engine 
torque, calculated by means of the torque model, even with 

3 0 discharge valves having highly advanced opening. This makes it 

possible to improve, especially simplify, the application of 
the torque model. 

Advantageous further developments and improvements of the 
35 method indicated in the main claim are rendered possible by 
the measures specified in the dependent claims. 
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It is advantageous, furthermore, if the efficiency is 
additionally determined as a function of the charge. Taking 
the charge into account improves the precision of the torque 
model even further- 

5 

An additional advantage consists in the fact that the 
deviation between an optimum ignition angle and an 
instantaneous ignition angle is selected as the variable 
characterizing the combustion center point. When the model is 
10 inverted, it allows both the calculation of a setpoint 

ignition angle at a given setpoint torque and given charge and 
the calculation of a setpoint charge at a given setpoint 
torque and given basic ignition-angle efficiency, and also the 
calculation of the actual torque. 

15 

A further advantage is that the efficiency is divided into a 
first partial efficiency and a second partial efficiency, the 
first partial efficiency being determined as a function of the 
variable characterizing the combustion center point and the 

2 0 second partial efficiency as a function of the variable 

characterizing the opening instant of the discharge -side 
gas-exchange valve. In this way, a simplified realization is 
possible if no extremely advanced opening of the 
discharge- side gas-exchange valve occurs. The efficiency may 
25 then be determined as product of the two partial efficiencies. 

Further advantages come to light from the following 
description of the exemplary embodiment. 

3 0 Brief Description of the Drawing 

Hereinafter, the present invention is explained in greater 
detail on the basis of the specific embodiments depicted in 
the drawing. The Figures show: Figure 1 a flow chart of a 
35 first specific embodiment of the utilized model; Figure 2 a 
general diagram of an engine control in which the sketched 
model is used; Figure 3 a further exemplary embodiment of the 
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torque model according to the present invention for 
determining an actual torque. Figure 4, by way of example, a 
map for determining the efficiency of the conversion of the 
chemical into mechanical energy; and Figure 5 an inverse 
5 torque model for calculating a setpoint ignition angle. 

Description of the Exemplary Embodiments 

Figure 1 shows a flow chart of a first exemplary embodiment 
10 for improving the torque model. Figure 1 describes a program 
of a microcomputer, the individual elements of the 
representation in Figure 1 representing programs, program 
steps or program parts, while the arrows describe the flow of 
information. Shown is the calculation of the actual basic 
15 torque, i.e., the torque that comes about when setting the 

basic ignition angle, which is selected from a characteristics 
map as a function of engine speed and load. 

It is essential that in the model shown in Figure 1 a 
2 0 correction of the optimum torque value takes place in which 

the inert -gas rate and the charge movement are considered and 
the optimum ignition angle value is corrected taking the 
inert-gas rate and the charge movement into account, [and] 
thus the instantaneous working point. 

25 

In a first ignition map 100, a value mioptlln for the optimum 
torque is read out as a function of engine speed nmot and 
actual charge rl, which is determined from the measured 
air-mass flow taking a suction manifold model into account. 

30 The ignition-map values of ignition map 100 are determined 

under defined, optimal marginal conditions, in particular when 
the lambda value has a standard value (1, for example) , an 
optimum ignition angle is set and efficiency etarri for the 
valve overlap (inert -gas rate and charge movement) has a 

35 standard value (for example, 1) . The optimum torque value is 
multiplied in a multiplication point 102 with efficiency 
etarri, which describes the deviation with respect to the 
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valve overlap from the standard value. Efficiency value 
etrarri is formed in ignition map 104 as a function of signals 
representing an inert -gas rate by internal and external 
exhaust-gas recirculation, and the charge movement- Proven 
5 useful has a signal rri for the internal and external 

inert-gas rate, which is calculated as a function of the 
setting of the exhaust -gas return valve and the settings of 
the intake and discharge valves. The inert-gas rate describes 
the portion of the inert gas of the total aspirated gas mass . 

10 Another way of calculating the inert -gas rate is based on the 
temperature of the recirculated exhaust-gas flow, lambda, the 
instantaneous air charge and the exhaust-gas pressure. A 
signal wnw, which represents the opening angle (relative to 
the crankshaft or the camshaft) of the intake valve, has shown 

15 to be suited for taking the charge movement into account . In 

other exemplary embodiments, the position of a charge -movement 
flap or a variable representing the lift and the opening phase 
of the intake valves is utilized. 

2 0 As a function of these instantaneous variables, efficiency 

etarri is determined, which describes the deviations in the 
torque value, attributable to the inert gas and the charge 
movement, from the torque value ascertained under standard 
conditions, on the basis of which ignition map 100 was 
25 determined. Optimum torque value mioptll, formed by the 

correction in multiplication point 102, is multiplied in a 
further multiplication point 106 by lambda efficiency etalam. 
It is ascertained in a characteristic curve 108 as a function 
of the instantaneous exhaust-gas composition lambda. The 

3 0 result is an optimum torque value miopt, which takes the 

instantaneous operating state of the internal combustion 
engine into account as well as its deviation from the standard 
values, which is used in determining the optimum torque 
values. Therefore, miopt is the optimum value of the indicated 
35 torque at an optimum ignition angle. To form basic torque 
mibas, from which the instantaneous torque may then be 
derived, the basic ignition angle setting relative to the 
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optimum ignition angle setting must therefore be taken into 
consideration. This occurs in multiplication point 110 where 
optimum torque value miopt is corrected by ignition-angle 
efficiency etadzw. 

5 

Ignition-angle efficiency etadzw is formed from a 
characteristic curve 112 as a function of the deviation, 
formed in 114, between basic ignition angle zwbas and optimum 
ignition angle zwopt . Efficiency etadzw thus determines the 

10 effect the deviation of the basic ignition angle from the 

optimum ignition- angle value has on the torque of the internal 
combustion engine. The basic ignition angle corresponds to the 
ignition angle that is read out from a basic ignition angle 
ignition map as a function of engine speed and load. It does 

15 not necessarily correspond to the actually set ignition angle 
that must be taken into account, in the form of an additional 
efficiency, when determining the actual torque as a function 
of the basic torque. The optimum ignition-angle value is read 
out from ignition map 116 as a function of engine speed nmot 

20 and charge rl. The read-out optimum ignition- angle value 
zwoptlln is conducted to a node 118 where this value is 
corrected by correction value dzworri. This correction value 
is generated in an ignition map 120 as a function of inert -gas 
rate rri, signal wnw for taking the charge movement into 

25 account, and of the signals, engine speed nmot and 

instantaneous charge rl representing the instantaneous 
operating point* Optimum ignition- angle value zwoptll 
corrected in this manner is corrected in an additional node 
122 by a correction value dzwola. It represents a lambda 

30 correction value, which is ascertained in ignition map 124 as 
a function of lambda and instantaneous charge RL. The 
corrections in nodes 118 and 122 are preferably implemented by 
addition. This approach has proven suitable since the 
instantaneous operating point of the internal combustion 

35 engine must then be considered in the correction values as 

well. The value zwopt resulting in 122 represents the optimum 
ignition- angle value which is compared in node 114 to the 
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basic ignition angle. 

The described model satisfactorily controls also the 
aforementioned operating states having high inert-gas rates 
5 and charge movements and small to medium fresh-air charges- It 
is essential in this context that a correction of the optimum 
ignition angle and a correction of the optimum torque take 
place in the torque model. The optimum ignition angle is made 
up of a basic value as a function of load and engine speed, 

10 and additive corrections as a function of the inert -gas rate, 
the settings of the intake and discharge valves or the valve 
overlap, of lambda, load and engine speed. The optimum torque 
is made up of a basic variable as a function of load and 
engine speed, and multiplicative corrections as a function of 

15 the inert -gas rate and the settings of the intake and 

discharge valves, of lambda and the ignition- angle efficiency 
with respect to the basic ignition angle. 

To determine the model parameters, a software tool is used, 
20 which is able to optimize the model parameters from the 

measured input variables and the measured torque, in such a 
way that the error square remains as small as possible across 
the entire measuring points. An example of such a software 
tool is known from DE 197 45 682 Al . 

25 

Basic torque mibas ascertained by means of the model is 
processed further in a variety of ways. Taking the efficiency 
of the actual ignition-angle setting into account, the 
instantaneous torque is calculated. Another evaluation 
30 consists in determining the ignition-angle setting, the 

difference between setpoint torque and basic torque being able 
to be utilized to correct the ignition-angle setting. 

The model described in Figure 1 shows the calculation of the 
35 instantaneous torque from various performance quantities. By 
reversing the model, the model, analogously to the model of 
the related art mentioned in the introduction, is also used to 
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determine the actuating variables (such as ignition angle, 
lambda, etc.), as a function of the setpoint torque value or 
of the deviation between setpoint torque and basic torque or 
instantaneous torque . 

5 

The model may also be formulated as follows; 

mibas = f 1 (nmot , rl) *f 21 (rri) *f 22 (wnw) *f 3 (lambda) *f 4 (zwopt- 
zwbas) 

10 

or, related to the instantaneous torque: 

miist = f 1 (nmot, rl) *f21 (rri) *f 22 (wnw) *f3 (0) *f4 (zwopt-zwist) 

15 By reversing the model, it is then possible to derive 

actuating variables such as an ignition-angle setpoint value 
zwsoll : 

zwsoll = zwopt - f4-^[misoll/(f21(rri)*f22(wnw)*f3(0) )] 

20 

The ignition maps and characteristic curves used to calculate 
the model are determined for each engine type within the 
framework of the application, possibly using the 
afore -mentioned software tool. 

25 

Figure 2 shows a cJontrol unit 400, which includes an input 
circuit 402, an output circuit 404 and a microcomputer 4 06. 
These components are connected to a bus system 408. The 
performance quantities to be evaluated for the engine control, 

30 which are detected by measuring devices 418, 42 0 through 424, 
are conveyed via input lines 410 and 412 through 416. The 
performance quantities required to calculate the model are 
shown above. The measured and possibly processed 
performance -quantity signals are then read in by the 

35 microcomputer via bus system 408. In microcomputer 406 itself, 
and there in its memory, the commands used for the model 
calculation are stored as computer program. This is symbolized 
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in Figure 2 by 426. The model results, which, if appropriate, 
are processed further in other programs (not shown) , are then 
conveyed from the microcomputer to output circuit 4 04 via bus 
system 408- Output circuit 404 thereupon outputs trigger 
5 signals as actuating variables for setting the ignition angle 
and the air supply, for example, and measured variables, such 
as instantaneous torque Mligf 

Figure 3 shows a flow chart for a torque model, using the 
10 aforementioned correlations for a basic operating point. In 
addition to the difference from optimum ignition angle and 
basic ignition angle (instantaneous ignition angle) , the 
variable of adjustment angle wnwao of the camshaft is also 
entered in block 510 for calculating an extended 
15 ignition-angle efficiency, this variable characterizing the 
opening instant of the discharge- side gas-exchange valve. 

As mentioned previously, an efficiency is generated in map 500 
as a function of engine speed Nmot and instantaneous charge 

20 rl, and in characteristic curve 502 as a function of inert-gas 
rate rri, the efficiency being multiplied in multiplication 
point 5 04 by the optimum torque value, which was formed in map 
500 and is optimal under standard conditions. In a 
multiplication point 506, optimum torque value mioptll 

25 determined therefrom is multiplied by a lambda efficiency, 

which is generated in characteristic curve 508 as a function 
of basic lambda value lambas, to be set in the given operating 
point without corrections possibly specified from the outside. 
The multiplication result in 506 is an optimum torque value 

30 miopt and is subjected to a further multiplication in 508 in 
which optimum torque value miopt is multiplied with 
ignition-angle efficiency etazwist formed in 510. The result 
is basic torque mibas for the instantaneous operating point. 
Analogously to the afore-mentioned representation, the 

35 ignition-angle efficiency is determined as a function of the 
difference, formed in 512, between optimum ignition angle 
zwopt and basic ignition angle zwbas as well as the directly 
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supplied adjustment angle wnwa of the camshaft. The optimum 
ignition-angle value is formed in a map 514 as a function of 
engine speed and charge; this optimum ignition angle is 
corrected in a summing point 516 as a function of a correction 
5 value determined as a function of charge movement LB, 

inert -gas rate rri, engine speed Nmot and charge rl. In a 
further correction point 518, the corrected optimum ignition 
angle is corrected by a correction value that is dependent on 
a lambda value and which is generated in characteristic curve 
10 52 0 as a function of the basic lambda value. Optimum 

ignition-angle value zwopt corrected in this manner is 
evaluated so as to generate the ignition-angle efficiency in 
510 and to determine the instantaneous torque. 

15 The torque model realized according to the flow chart in 

Figure 3 is thus also suitable for internal combustion engines 
having discharge -side gas-exchange valves or discharge valves 
that have advanced or very advanced opening* On the basis of 
basic torque mibas, it provides the instantaneous torque as 

20 indicated engine torque with high accuracy even with discharge 
valves having advanced or very advanced opening and retarded 
ignition angles . An advanced opening of a discharge valve may 
be achieved, for example, by a corresponding phase shift of 
the discharge camshaft, 

25 

A discharge valve is considered to have advanced opening in 
this context if it opens prior to reaching bottom dead center 
of ignition- phase ignition-BDC, and to have very advanced 
opening if it opens even less than a crank angle, or less than 
3 0 a 120** phase shift of the discharge camshaft after TDC of 
ignition phase ignition-TDC, and thus more than a 60° crank 
angle prior to reaching BDC of ignition phase ignition-BDC. 

Decisive for the conversion of chemical into mechanical energy 
35 and thus for generating the indicated torque is the 

utilization of the heat energy that is released in the course 
of combustion. 
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Disregarding the wall-heat losses in the combustion chamber, 
the opening instant of the individual discharge valve is also 
decisive for this utilization of the heat energy, in addition 
to the so-called combustion center point and the progression 
5 of the combustion, which is known as fibe function. 

The combustion center point is defined in this context as the 
crank angle at which half of the combustion energy has been 
converted. The fibe function describes the temporal 
10 conversion of the chemical energy over the crank angle. 

A very advanced opening of the individual discharge valve 
leads to a rapid drop in the pressure in the combustion 
chamber, at a time when the piston is still noticeably moving 
15 downward, and thus leads to a lower averaged indicated torque. 
The portion of the "lost" torque not only depend on the 
opening instant of the individual discharge valve, that is to 
say, the crank angle at which the individual discharge valve 
opens, for instance, but also on other variables. 

20 

The combustion center point, in particular, plays a decisive 
role. In a retarded combustion center point, which is 
characterized by retarded ignition angles, for instance, the 
pressure is much higher at a crank angle of 12 0 degrees after 
25 top dead center of ignition phase ignition-TDC, for instance. 

The working share of the mechanical energy provided during the 
second half of the downward travel of the piston is greater 
with retarded ignition angles . 

30 The smaller the charge in the combustion chamber, the lower 
also the combustion peak pressure and the faster the 
combustion- chamber pressure drop during the downward travel of 
the piston. With very small loads, the combustion- chamber 
pressure falls below the ambient pressure, so that in this 

35 case an advanced opening of the individual discharge valve 
leads to an increase in the indicated torque, since the 
negative portion of the expansion work is left out . 
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However, with retarded ignition angles, a maximum indicated 
torque is obtained even with small charges upon opening of the 
individual discharge valve at bottom dead center of ignition 
phase ignition-BDC or shortly before bottom dead center of 
5 ignition phase ignition-BDC . However, simulation results have 
shown that, with respect to the relative variables, there is 
hardly any change as a function of charge rl . Therefore, 
charge rl is not necessarily required for considering the 
losses of an individual discharge valve having advanced 
10 opening. However, its consideration would improve the accuracy 
of the torque model slightly more. For this reason, Figure 3 
provides charge rl as additional input variable in block 510, 
as a dashed line. 

15 The structure of the torque model described in the not 
pre-published German patent application 101 49 477.7 is 
broadened here according to the specific embodiment of Figure 
3. In this context, the currently defined ignition-angle 
efficiency etadzw = f (dzw) is expanded in its importance. It 

20 is given an additional dependency on crank angle wnwao at 
which the individual discharge valve opens. It therefore 
describes the overall efficiency of the conversion of chemical 
into mechanical energy, taking also those losses into account 
that result from an advanced opening of the individual 

25 discharge valve. 

The flow chart according to Figure 3 shows an exemplary 
realization. In contrast to what is known, a characteristic 
map KFETADZW is now realized in block 510, which determines 

3 0 the overall efficiency etadzw of the conversion of chemical 
energy into mechanical energy, for instance from crank angle 
wnwao at which the individual discharge valve opens und which 
corresponds to the adjustment angle of the discharge camshaft, 
and deviation dzw of the ingnition angle from its optimum 

35 value zwopt at which the greatest induced torque comes about. 

Characteristic map KFETADZW is shown in Figure 4 as result of 
a simulation. The reduction of the indicated torque, i.e., the 
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reduction of overall efficiency etwadzw, is [word/s missing^] 
with simultaneous retarding of the ignition, that is, for 
increasing values of deviation dzw of the ignition angle from 
its optimum value zwopt, and advancing of the opening instant 
5 of the individual discharge valve, i-e,, reduction of crank 
angle wnwo at which the individual discharge valve opens, 
relative to top dead center of ignition phase igriition-TDC* 

Instead of the deviation of the ignition angle from its 
10 optimum value zwopt, it is also quite possible to utilize the 
combustion center point or some other variable as input 
variable of characteristic map KFETADZW, which describes the 
position of the combustion over the crank angle* 

15 Instead of adjustment angle wnwao of the discharge camshaft, 

it is also possible to use as input variable of characteristic 
map KFETADZW some other variable characterizing the opening 
instant of the individual discharge valve. 

20 In a simplified realization, in the event that no extremely or 
very advanced opening of the discharge valve occurs, the 
possibility exists, at the expense of accuracy, to reduce 
characteristic map KFETADZW into the product of two 
characteristic lines. This corresponds to dividing the overall 

25 efficiency into two partial efficiencies. A first partial 
efficiency is determined as a function of a variable 
characterizing the combustion center point, i.e., the 
deviation of the ignition angle from its optimum value zwopt, 
for example. A second partial efficiency is determined as a 

30 function of the variable characterizing the opening instant of 
the individual discharge valve. The first partial efficiency 
in these examples therefore is the current pure ignition- angle 
efficiency, and the second partial efficiency in this example 
is the efficiency of the opening instant of the respective 

^ Translator's Footnote? A word or sentence fragment is missing in the 
German source text . 
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discharge valve. 



For the sake of simplification, the influence of charge rl on 
the overall efficiency will not be considered in this 
5 examination. 

Taking charge rl into account would result in a 
three-dimensional characteristic map. Charge rl may be a 
fresh-air charge or an exhaust -gas -enriched overall charge 

10 provided an external exhaust -gas recirculation is present or 
the internal exhaust-gas recirculation is of importance. In 
this case, in addition to charge rl, inert -gas rate rri could 
also be utilized as input variable for characteristic map 
KFETADZW, which would then even become four-dimensional or 

15 five- dimensional . To be able to consider all influences on the 
overall efficiency, charge movement lb could additionally be 
included as well. In Figure 3, inert -gas rate rri and charge 
movement lb, since optional, are likewise represented as input 
variables of characteristic map 510 in the form of dashed 

20 lines. 

Using the improved torque model according to the flow chart of 
Figure 3 achieves high precision in the calculation of the 
indicated torque even in individual discharge valves that have 
25 very advanced opening and retarded ignition angles. An 

advanced opening of the individual discharge valve may be 
desired to push the hot combustion gases into the exhaust-gas 
system early on, thereby substantially accelerating the 
heating of the catalytic converter. 

30 

The aforementioned torque model according to Figure 3 must be 
inverted to calculate the ignition setpoint ignition angle. 
Figure 5 shows a corresponding block diagram. 

35 In a block 1000, setpoint value etazwsoll for the expanded 
ignition- angle efficiency or the overall efficiency is 
determined by forming the quotient from setpoint torque misoll 
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to be set and optimum torque miopt, using the following 
formula: 

etazwsoll = misoll/miopt . 

5 

Inverse characteristic map 1100, also called KFDZWETA, then 
receiveses as input variables setpoint value etazwsoll for the 
overall efficiency and crank angle wnwao at which the 
individual discharge valve opens. Once again, it is optionally 

10 possible to transmit to characteristic map KFDZWETA as input 

variable a value for charge rl, possibly taking inert -gas rate 
rri and charge movement lb into account. A deviation dzwsoll 
of setpoint ignition angle zwsoll from optimum ignition- angle 
value zwop will then result as output variable. In order to 

15 determine setpoint -ignt ion angle zwsoll, deviation dzwsoll of 
the setpoint -ignition angle is then deducted from optimum 
ignition-angle value zwopt in a subtraction point 1200. 
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